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Abstract 

Sampled-data  analysis  of  converters  has  been  a  topic  of  investigation  for  the  past  two  decades. 
However,  this  powerful  tool  is  not  widely  used  in  control  loop  design  or  in  closed-loop  perfor¬ 
mance  validation.  Instead,  averaged  models  are  typically  used  for  control  loop  design,  while 
detailed  simulations  are  used  for  validating  closed-loop  performance.  This  paper  makes  several 
contributions  to  the  sampled-data  modeling  and  analysis  of  closed-loop  PWM  DC-DC  convert¬ 
ers,  with  the  aim  of  increasing  appreciation  and  use  of  the  method.  General  models  are  presented 
in  a  unified  and  simple  manner,  while  removing  simplifying  approximations  present  in  previous 
work.  These  models  apply  both  for  current  mode  control  and  voltage  mode  control.  The  gen¬ 
eral  models  are  nonlinear.  They  are  used  to  obtain  analytical  linearized  models,  which  are  in 
turn  employed  to  obtain  local  stability  results.  Detailed  examples  illustrate  the  modeling  and 
analysis  in  the  paper,  and  point  to  situations  in  which  the  sampled-data  approach  gives  results 
superior  to  alternate  methods.  For  instance,  it  is  shown  that  the  sampled-data  approach  will 
reliably  predict  the  (local)  stability  of  a  converter  for  which  averaging  or  simulation  predicts 
instability. 


1  Introduction 

Sampled-data  analysis  of  converters  has  been  a  topic  of  investigation  for  the  past  two  decades 
[1,  2,  3].  Sampled-data  models  of  converters  are  accepted  to  be  more  accurate  than  their  averaged 
counterparts.  This  can  be  important  especially  in  studying  dynamical  behaviors  that  averaging 
is  known  not  to  captnre  well,  such  as  subharmonic  instability  [4],  chaotic  phenomena  [5,  6,  7,  8], 
and  steady-state  DC  offset  [9].  However,  the  powerful  sampled-data  approach  is  not  widely  used 
in  control  loop  design  or  in  closed-loop  performance  validation.  Instead,  averaged  models  are 
typically  used  for  control  loop  design,  while  detailed  simulations  are  used  for  validating  closed-loop 
performance.  This  paper  makes  several  contributions  to  the  sampled-data  modeling  and  analysis 
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of  closed-loop  PWM  DC-DC  converters,  with  the  aim  of  increasing  appreciation  and  use  of  the 
method.  General  models  are  presented  in  a  unified  and  simple  manner,  while  removing  simplifying 
approximations  present  in  previous  work. 

Because  the  dynamical  behavior  differs  for  continuous  conduction  mode  (CCM)  and  discontin¬ 
uous  conduction  mode  (DCM),  the  general  model  in  CCM  differs  from  that  in  DCM.  However,  in 
either  case  the  model  applies  both  for  current  mode  control  and  voltage  mode  control.  Because  of 
this,  the  modeling  procedure  presented  here  is  considered  to  unify  previous  disparate  contributions 
for  current  mode  and  voltage  mode  control. 

The  general  models  are  nonlinear.  They  are  used  to  obtain  analytical  linearized  models,  which 
are  in  turn  employed  to  obtain  local  stability  results.  Detailed  examples  illustrate  the  modeling 
and  analysis  in  the  paper,  and  point  to  situations  in  which  the  sampled-data  approach  gives  results 
superior  to  alternate  methods.  For  instance,  it  is  shown  that  the  sampled-data  approach  will  reliably 
predict  the  (local)  stability  of  a  converter  for  which  averaging  or  simulation  predicts  instability. 

It  is  important  to  clarify  at  the  outset  the  differences  between  the  contributions  to  modeling 
in  this  paper  in  relation  to  past  contributions.  In  the  next  paragraph,  the  relation  between  the 
present  work  and  the  original  work  on  sampled-data  modeling  and  analysis  of  closed-loop  PWM 
DC-DC  converter  operation  is  summarized.  This  is  followed  by  a  summary  of  the  relation  between 
this  work  and  the  literature  on  sampled-data  modeling  and  analysis  of  the  power  stage  of  a  PWM 
DC-DC  converter. 

Lee,  Iwens,  Yu  and  Triner  [1]  and  Verghese,  Elbuluk  and  Kassakian  [2]  are  two  seminal  con¬ 
tributions  to  sampled-data  modeling  and  analysis  of  DC-DC  converters.  These  authors  developed 
general  models  that  apply  to  closed-loop  converter  operation.  Thus,  their  results  apply  also  to  the 
modeling  of  the  power  stage  of  a  converter.  In  relation  to  these  papers,  the  present  contribution 
gives  more  detail,  unifies  the  modeling  of  current  mode  control  and  voltage  mode  control,  and  gives 
an  analytical  linearized  sampled-data  model  (not  given  explicitly  in  the  indicated  references).  In 
addition,  the  models  given  here  incorporate  the  exact  system  trajectory  between  switching  instants. 

Tymerski  [10,  11]  developed  general  detailed  sampled-data  models  of  the  power  stage  of  a  PWM 
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DC-DC  converter.  He  also  obtained  an  analytical  linearized  sampled-data  model.  These  papers  are 
similar  to  the  present  work,  especially  in  that  they  reflect  the  exact  trajectory  between  switching 
instants  and  inclnde  an  analytical  linearized  sampled-data  model.  The  main  difference  between  the 
work  in  [10,  11]  and  the  present  paper  is  that  here  the  closed-loop  converter  is  considered,  while  in 
[10,  11]  the  power  stage  is  considered. 

The  paper  proceeds  as  follows.  In  Section  2,  a  general  block  diagram  model  valid  in  continnons 
condnction  mode  is  given  and  nsed  to  obtain  a  nonlinear  sampled-data  model  and  an  associated 
linearized  model.  The  linearized  model  is  used  to  stndy  orbital  stability,  and  to  derive  the  andio- 
susceptibility  and  outpnt  impedance.  In  Section  3,  analogous  models  are  given  for  discontinnons 
conduction  mode.  In  Section  4,  four  examples  are  given  that  illustrate  application  of  the  models  of 
the  paper  to  sample  converter  circuits.  In  addition,  these  examples  demonstrate  advantages  of  the 
sampled-data  analysis  over  averaging  and  detailed  simulation  in  validating  closed-loop  performance. 
Conclnsions  are  collected  in  Section  5. 

2  Continuous  Conduction  Mode  (CCM) 

2.1  Block  Diagram  Model 

In  this  subsection,  a  block  diagram  model  for  the  PWM  converter  in  CCM  is  proposed.  This  model 
is  shown  to  apply  to  cnrrent  mode  control  and  voltage  mode  control. 

Consider,  for  example,  a  buck  converter  under  cnrrent  mode  control  and  voltage  mode  control, 
shown  in  Fig.  1  and  Fig.  2,  respectively.  They  have  different  control  schemes.  However,  they  share 
the  same  characteristics:  switching  instants  are  determined  by  how  a  ramp  signal  h{t)  intersects 
with  a  feedback  signal  y{t).  This  motivates  a  nnifled  model  for  both  current  mode  control  and 
voltage  mode  control. 

The  proposed  block  diagram  model  for  the  PWM  converter  in  CCM  is  shown  in  Fig.  3.  In  the 
diagram,  Ai,A2  G  Bi,B2  G  C,Ei,E2  G  R^^^,  and  D  G  R  are  constant  matrices, 

X  G  R^,  y  =  Cx  Du  G  R  are  the  state  (of  power  stage  and  controller)  and  the  feedback  signal, 
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Figure  1:  Buck  converter  under  cnrrent  mode  control 


Figure  2:  Buck  converter  under  voltage  mode  control 
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respectively.  The  source  voltage  is  Vg]  the  output  voltage  is  Vq-  The  notation  Vr  denotes  the 
reference  signal,  which  could  be  a  voltage  or  current  reference.  The  reference  signal  Vr  is  allowed 
to  be  time-varying,  although  it  is  constant  in  most  applications.  The  signal  h{t)  is  a  T-periodic 
ramp  with  h{0)  =  Vi  and  h{T~)  =  V^.  The  clock  has  the  same  freqnency  fg  =  1/T  as  the  ramp. 
This  frequency  is  called  the  switching  frequency.  The  two  stages  in  each  clock  period  in  CCM  are 
denoted  by  Si  and  52.  The  system  is  in  5i  immediatedly  following  a  clock  pulse,  and  switches  to 
S2  at  instants  y{t)  =  h{t). 

h{t)  =  Vi  +  {Vh-Vi){^xnod  1) 
clock 

_  y  =  Cx  +  Du 


u={  ) 

Vr 


Switch  to  Si  or  S2 


\  X 
[  Vo 

=  Aix  -|-  Biu 
=  Eix 

52  :  1 

\  X 

=  A2X  -1-  B2U 

[  Vo 

=  E2X 

Switching 

Decision 


Figure  3:  A  block  diagram  model  for  PWM  converters  in  CCM 

Assnme  the  switch  and  the  diode  in  the  PWM  converter  are  ideal,  so  that  there  is  no  voltage 
drop  when  they  are  on.  Most  PWM  converters  can  be  modeled  by  Fig.  3.  The  PWM  converter 
under  current  mode  control  fits  this  model  exactly,  with  h{t)  denoting  the  slope-compensating 
ramp.  For  operation  in  voltage  mode  control,  the  system  is  switched  between  stages  when  the 
ramp  signal  h{t)  intersects  with  the  feedback  signal  y.  One  of  the  switchings  generally  has  the 
same  frequency  as  the  clock.  Therefore  the  model  in  Fig.  3  is  also  good  for  the  PWM  converter 
under  voltage  mode  control.  Many  other  control  schemes  (e.g.,  average  current  mode  control)  also 
fit  the  model  in  Fig.  3.  Typical  waveforms  in  cnrrent  and  voltage  mode  control  are  shown  in  Fig.  4 
and  Fig.  5,  respectively.  Note  that  in  Fig.  4,  the  ramp  has  positive  slope,  instead  of  negative  slope 
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as  commonly  seen  in  most  literatures,  in  order  to  be  consistent  with  the  case  of  voltage  mode 
control. 


S-)  S2  Si  Sg  Si  S2 


Figure  4:  Waveforms  of  a  PWM  converter  in  CCM  under  current  mode  control 


h(t) 

y(t) 


Switch 

Si  S2  Si  S2  Si  S2 


Figure  5:  Waveforms  of  a  PWM  converter  in  CCM  under  voltage  mode  control 


2.2  Nonlinear  Sampled-Data  Model 

The  switching  action  is  essentially  discrete.  To  arrive  at  a  sampled-data  model,  the  operation  of 
the  PWM  converter  within  the  cycle,  t  €  [nT,  (n  +  1)T),  is  considered.  Generally  in  the  PWM 
converter,  the  switching  frequency  is  sufficiently  high  that  the  variations  in  Vg  and  Vr  in  a  cycle 
can  be  neglected.  Thus,  take  u  =  {vs,Vr)'  G  to  be  constant  within  the  cycle,  and  denote  its 

value  by  Un  =  {vsniVm)' ■  (The  notation  Vsn-,  instead  of  Vg^n-,  is  used  for  brevity.  This  applies  to 
other  variables.)  Let  Xn  =  x{nT)  and  Von  =  Vo{nT).  The  sampled-data  dynamics  which  maps  Xn 
to  Xn+i  is  derived  next. 

Denote  by  nT  +  dn  the  switching  instant  within  the  cycle  when  y{t)  and  h{t)  intersect.  Then 
y{nT  -I-  dn)  =  h{nT  -|-  dn),  and 

Si  :  I  ~  foi  t£[nT,nT  +  dn)  (1) 
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52 


X 


A2X  +  B2U 
E2X 


for  t  G  [nT  +  dn,  {n  +  l)r) 


(2) 


O 


The  two  matrices  Ei  and  E2  need  not  be  the  same.  For  example,  they  can  differ  if  the  equivalent 
series  resistance  (ESR)  Rc  7^  0.  When  they  differ,  the  ontput  voltage  is  discontinuous.  An  example 
of  a  discontinuous  output  voltage  waveform  is  shown  in  Fig.  6.  In  most  applications,  the  output 
voltage  of  interest  is  the  maximum,  minimum,  or  average  voltage.  So  in  the  following,  E  is  used  to 
denote  either  Ei,  E2,  or  {Ei  +  i?2)/2. 

a> 
re 

"o 
> 

3 
a 
3 
O 

Ns 


nT  (n+1)T  t 


Ei  x(nT) 


E2X(nT) 


Figure  6:  A  discontinuous  output  voltage  waveform 


Using  the  dynamics  (1)  and  (2)  and  the  assumption  that  u  =  {vs,Vr)  is  constant  within  the 
cycle,  one  readily  obtains  the  sampled-data  dynamics  (3)-(5)  of  the  PWM  converter  in  CCM. 
The  dynamical  equation  (3)  below  is  augmented  with  the  constraint  (4),  corresponding  to  the 
intersection  condition  y{nT  +  dn)  =  h{nT  +  dn)  defining  dn- 


Xn+l  —  fiXnjUmdn) 

fdn  pT 

=  +  /  e^^^^-^UaB2Un  (3) 

-/  0  J  dn 

g{xn,Un,dn)  =  y{nT  +  dn)  -  h{nT  +  dn) 


=  '^^daBiUn)  +  Dun  -  h{dn) 

Jo 


=  0 

Von  —  EXn 


(4) 

(5) 
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2.3  Steady-State  Analysis 


In  the  PWM  converter,  the  steady-state  operating  condition  is  a  periodic  solution,  not  an  equi¬ 
librium  point  as  depicted  in  the  averaging  method.  A  T-periodic  solution  for  system  in  Fig.  3 
corresponds  to  a  fixed  point  for  the  sampled-data  model.  A  periodic  solution  x^{t)  is  sampled  as 
x°(0).  Let  the  fixed  point  in  the  sampled-data  dynamics  (3)-(5)  be  {xn,Un,dn)  =  {x^{0),u,d), 
where  u  =  (I4,  Vr)'-  Then  this  fixed  point  satisfies 

=  fix^iO),u,d)  (6) 

5(x°(0),«,d)  =  0  (7) 


where  the  functions  /(•)  and  g{-)  are  given  in  Eqs.  (3)  and  (4).  These  N  +  1  nonlinear  equations 
(Eqs.  (6)  and  (7))  in  A  -|-  1  unknowns  (x°(0)  and  d)  can  be  solved  by  Newton’s  method  [1].  After 
obtaining  x°(0)  and  d,  a  periodic  solution  x^{t)  is  obtained: 


x«(t) 


gAit^0('Q)  _|_  dcrBiu  for  t  G  [0,  d) 

<  _j_  fortG[d,  T) 

x^{t  mod  T)  for  t  >  r 


(8) 


A  typical  periodic  solution  x^{t)  is  shown  Eig.  7.  Erom  Eqs.  (1)  and  (2),  x^{d  )  =  Aix^{d)  +  Biu 
and  x^{d~^)  =  A2X^{d)  +  B2U. 


Eigure  7:  A  typical  periodic  solution  x^{t)  of  a  PWM  converter  in  state  space 


2.4  Existence  of  T-Periodic  Solutions 

Assume  all  of  the  eigenvalues  of  Ai  and  A2  are  in  the  open  left  half  of  the  complex  plane,  then 
the  matrices  /  —  Q^2(T-d) ^Aid  j  _  ^Aid^A2(T-d)  invertible  [12].  Prom  Eq.  (6),  x°(0)  can  be 


8 


expressed  as 


x°(0)  =  (/_e^2(T-d)gAid^-i^gA2(T-d)  e^^^'^-^UaB2u)  (9) 

Jo  Jd 

Similarly,  x^{d)  can  be  expressed  as  a  function  of  d,  denoted  as  Xs{d): 


Xsid)  =  {I r  e^^^'^-^UaB2U+  e^^^‘^-^UaBiu)  (10) 

Jd  Jo 

The  function  Xs(-)  has  the  following  property: 

X,(0)  =  -A:^^B2u  (11) 

XsiT)  =  -Ai^Biu  (12) 

So  the  +  1  equations,  (6)  and  (7),  reduce  to  one  equation  in  one  unknown  d: 

y{d)  -  h{d)  =  CXs{d)  +  Du  -  h{d)  =  0  (13) 


The  next  result  gives  a  sufficient  condition  for  existence  of  a  fixed  point  in  the  sampled-data 
dynamics. 

Theorem  1  Assume  that  all  of  the  eigenvalues  of  A\  and  A2  are  in  the  open  left  half  of  the  complex 
plane.  If 

{CAf^B2U  -Du  +  h{{))){CAf^Biu  -Du  +  h{T^))  <  0  (14) 

then  there  exists  a  fixed  point  {xn,Un,dn)  =  (x^(0),  u,  d)  in  the  sampled-data  dynamics  (3)-(5). 
Proof:  From  Eq.  (13),  if 

(y(0)-h(0))(y(r-)-h(r-)) 

=  (CX,(0)  +Du-  h{0)){CXs{T)  +Du-  h{T-)) 

=  {CAf^B2U  -Du  +  h{0)){CAf^Biu  -Du  +  h{T-)) 

<  0 
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then  by  the  Intermediate  Value  Theorem,  there  exists  a  solution  d  €  [0,T)  satisfying  Eq.  (13). 
Hence  there  exists  a  fixed  point  {xn,Un,dn)  =  (x^(0),  u,  d)  in  the  sampled-data  dynamics  (3)-(5). 
□ 

Note:  Whenever  there  exists  a  T-periodic  solution  for  the  system  in  Fig.  3,  it  will  be  found  as 
a  fixed  point  in  the  sampled-data  model.  However,  a  fixed  point  in  the  sampled-data  model  may 
not  correspond  to  a  T-periodic  solution  in  Fig.  3.  To  explain  the  last  remark,  suppose  a  fixed  point 
{xn,Un,dn)  =  {x^{0),u,d)  is  obtained.  Then  a  solution  x^{t)  can  be  obtained  from  Eq.  (8).  An 
example  of  signals  y^{t)  :=  Cx^{t)  +  Du  and  h{t)  is  shown  in  Fig.  8.  The  solution  although 

periodic  with  y^{d)  =  h{d),  cannot  be  accepted  because  two  switchings  occur  during  t  €  (0,T). 


Figure  8:  Ineligible  signals  y^(t)  and  h(t)  in  voltage  mode  control 


2.5  Linearized  Sampled-Data  Dynamics 

The  sampled-data  dynamics  (3)-(5)  is  constrained  and  nonlinear.  Assuming  that  =  Cx^(d~)  — 

h(d)  7^  0,  the  Implicit  Function  Theorem  allows  one  to  linearize  (3)-(5)  at  the  fixed  point  (xn,  Un,  dn)  = 
{x^{0),u,d)  to  obtain  an  unconstrained  linear  sampled-data  dynamical  model.  Using  a  hat  ~  to 
denote  small  perturbations  (e.g.,  Xn  =  Xn  —  a:°(0)),  it  follows  that 

Xn+l  ^  ^Xn  +  TUn  ,  . 

V  -  Ft 

Uon  - 

where 

$  =  (  ^9  ^-1  dg 

dXn  ddn  Odn  dXn  (xn,Un,dn)=(x^{0),u,d) 
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=  e 


=  e 


A2(T-d)rj  _  ((^1  ~  M)x'^{d)  +  {Bi  -  B2)u)C X  Aid 
C{Aix^{d)+Biu)-h{d)  ^ 

M{T-d)(j-  _  (x°(d~)-a;°(d+))C 
^  CxO(d-)  -  h{d) 


m 


r  = 


df  ^9 

V  Q  J  / 


dUr,  ddr,  ddn  du. 


(Xn,Un,dn)  =  {xO{0),u,d) 

x^{d~)  —  x°(d+) 


"M  /  e^^^daBi  —  . 

Vo  C7xO(d-)  -  /i(d)  io 


pT—d 

{C  /  e^^^daBi  +  D))  +  / 

JO  JO 


Recall  that  Un  =  {vsn-iVm)' ■,  and  denote  F  =  [ri,r2].  Eq.  (15)  can  be  rewritten  as 


J^n+l  ~  “1“  +  r2'0rn 

Von  —  BXji 


(18) 


2.6  Stability  Analysis 

The  relevant  stability  notion  is  asymptotic  orbital  stability,  not  asymptotic  stability  of  an  equilib¬ 
rium  point  as  depicted  in  the  averaging  method.  In  the  power  electronics  literature,  the  PWM 
converter  is  generally  said  to  be  either  stable  or  unstable,  without  mentioning  orbital  stability  per 
se.  The  definition  of  asymptotic  orbital  stability  is  given  as  follows. 

Definition  1  (see,  e.g.,  [13])  Denote  by  7  the  elosed  orbit  generated  by  the  periodie  solution  x^{t). 
Then  x^{t)  is  asymptotically  orhitally  stable  if  there  is  a  6  such  that 

dist[a;(0),  7]  <  6  ^  ^lim  dist[x(t),7]  =  0 

where  dist[z,7]  is  defined  as  the  smallest  distance  between  the  point  z  and  any  point  on  7. 

The  orbital  stability  is  related  to  the  the  set  of  eigenvalues  of  ^>,  denoted  by  (t[$].  Recall  the 
following  basic  results  on  asymptotic  orbital  stability  of  the  periodic  solution  x^{t). 

Theorem  2  [13]  The  fixed  point  x^ (0)  of  the  system  (3)-(5)  is  asymptotically  stable,  or  equivalently 
the  periodic  solution  x^{t)  in  the  original  continuous-time  system  of  Fig.  3  is  asymptotieally  orhitally 
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stable  if  all  of  the  eigenvalues  of  $  are  inside  the  unit  cirele  of  the  complex  plane.  Moreover,  if 
x^{f)  is  asymptotically  orbitally  stable,  then  no  eigenvalues  of  ^  lies  outside  the  unit  circle. 

The  following  well  known  fact  from  linear  algebra  [14,  p.  667]  will  be  used  in  the  proofs  below. 

Lemma  1  If  A  and  B  are  conformable  reetangular  matriees,  then 

det[I  +  AB]  =  det[I  + BA] 


The  preceding  theorem  and  lemma  will  be  used  to  show  the  following. 


Theorem  3  If  the  periodic  solution  x^{t)  is  asymptotically  orbitally  stable,  then  the  following  in¬ 
equality  holds: 


Cx^{d+)  -  h{d) 
Cx^{d-)  -  h{d) 


^  gtr[A2-Ai](i-tr[A2]T 


(19) 


Equivalently,  if  (19)  does  not  hold,  then  the  periodic  solution  x^{t)  is  not  asymptotically  orbitally 
stable. 


Proof:  Snppose  the  periodic  solution  x^{t)  is  asymptotically  orbitally  stable.  Then  all  the  eigenval¬ 
ues  of  $  have  magnitude  less  than  or  equal  than  1.  Since  det[4>]  is  the  prodnct  of  the  eigenvalues 
of  4>,  it  follows  that 


|det[4>]|  = 


det[e^i''e^^(^-'')]  det[/  -  )-x°(d+))C 

^  ^  ^  CxO{d-)-h{d)  ' 


det[e^i''e^2(^-'')]  det[l  -  )-^°(^+))- 

Cx^{d-)  -  h{d)  ^ 


^  g-tr[A2-Ai]d+tr[A2]T 


Cx^{d+)  -  h{d) 


Cx^id-)  -  h{d) 


<  1 


where  the  second  line  of  the  calculation  follows  from  Lemma  1. 


□ 
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Remark:  Inequality  (19)  is  a  necessary  (but  not  sufficient)  condition  for  stability.  Generally  the 
switching  period  is  so  small  that  the  right  side  of  (19)  can  be  approximated  as  1,  resulting  in  a 
condition  that  resembles  a  well-known  stability  criterion  in  current  mode  control  [15,  for  example]: 


m2  —  me 
mi  —  me 


(20) 


where  mi  is  the  (positive)  slope  of  the  inductor  current  trajectory  during  the  on  stage  and  m2  is 
the  (negative)  slope  during  the  off  stage  using  a  linear  approximation  [16];  and  me  is  the  (negative) 
slope  of  the  compensating  ramp.  (Note  that  here  mi  >  0,  m2  <  0  and  me  <  0,  which  differs  from 
most  literatures’  notations:  mi  >  0,  m2  >  0  and  me  >  0.)  Theorem  3  differs  from  (20)  in  its  use 
of  instantaneous,  rather  than  approximate,  slope.  Also,  Cx^{t)  in  Theorem  3  is  composed  of  more 
state  signals  (up  to  N)  than  just  inductor  current. 

The  next  result  asserts  that  a  closed-loop  pole  is  zero  if  the  ramp  has  slope  CxP{d'^)  = 
C{A2x'^{d)  +  B2u). 

Theorem  4  If  the  ramp  slope  h{d)  =  Cx^{d'^),  then  the  matrix  has  a  zero  eigenvalue. 

Proof:  From  the  proof  in  Theorem  3,  it  can  be  seen  that  if  h{d)  =  Cx^{d~^),  then  det[$]  =  0  and 
hence  the  matrix  has  a  zero  eigenvalue.  □ 

Remark:  Generally,  a  ramp  with  slope  Cx^{d'^)  is  considered  to  produce  deadbeat  behavior 
(finite  settling  time)  [15,  p.  418].  This  requires  all  of  the  N  closed-loop  poles  to  be  at  0.  Theorem 
4  predicts  that  if  the  ramp  has  slope  Cx^{d'^),  then  one  eigenvalue  is  0.  Indeed,  examples  can  be 
constructed  showing  that  not  necessarily  all  the  eigenvalues  vanish.  Deadbeat  behavior  is  therefore 
not  guaranteed. 

Theorem  5  and  Gorollary  2  below  address  the  case  in  which  the  switching  frequency  goes  to 
infinity.  The  following  corollary  of  Lemma  1  is  used  in  the  proof  of  Theorem  5. 


Corollary  1  Let  a  and  b  be  column  vectors  of  the  dimension  N.  Then  one  eigenvalue  of  the  square 
matrix  ab'  is  b'a,  and  the  remaining  iV  —  1  eigenvalues  are  zero. 
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Proof: 


det[A/  —  ab'] 

=  A^det[/-^] 

A 

(21) 

=  det[l  — 

A 

(22) 

=  A^(l-^) 

A 

(23) 

=  A^-^(A-5'a), 

(24) 

where  the  second  line  follows  from  Lemma  1. 

□ 

Theorem  5  As  the  switching  frequency  fs  goes  to  infinity,  an  eigenvalue  of  $  asymptotically  ap¬ 
proaches  remaining  (N  —  1)  eigenvalues  approach  1. 


Proof:  As  the  switching  frequency  fs  goes  to  infinity,  both  and  e^^'^  approach  identity 


matrices.  Eq.  (16)  now  implies  that  $  approaches  (I  —  matrix 

(x^(d~)  —  x^(d~^))C  is  the  product  of  a  column  vector  and  a  row  vector,  so  that  Corollary  1  applies. 
Thus  this  matrix  has  an  eigenvalue  at  C(x^(d~) —  x^(d~^))  and  (A^—  1)  eigenvalues  at  0.  This  in  turn 


implies  that  the  matrix  I  —  ^  lias  an  eigenvalue  at  1  — 

Cx^{d~  )—h(d)  Cx^{d~  )—h(d)  Cx^{d~  )—h{d) 

and  (A”  —  1)  eigenvalues  at  1.  □ 

Remark:  Generally  the  averaged  model  is  reliable  for  sufficiently  high  switching  frequency  [9]. 
This  theorem  implies  that  the  stability  criterion  in  (20)  does  indeed  provide  a  necessary  condition 
for  stability  when  the  switching  frequency  is  sufficiently  high. 


Corollary  2  Assume  the  switching  frequency  goes  to  infinity  and  Cxf'{d  )  <  h{d)  (which  is  gen¬ 
erally  true).  If  (i)  Cx^{d'^)  <  Cx^{d~)  or  (ii)  h{d)  <  {d+)+x  {d  ))  ^  system  is  not 

asymptotically  stable. 


Proof:  Under  the  stated  assumptions,  it  can  be  proved  that  the  eigenvalue  greater 


than  1  under  condition  (i),  or  less  than  -1  under  condition  (ii).  Thus,  the  system  is  not  asymptot- 
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ically  stable.  □ 

Remark:  The  corollary  asserts  that  if  the  switching  frequency  is  high  and  if  the  slopes  of  the 
output  variable  ij  =  Cx  are  not  properly  set  in  the  design  (condition  (i))  or  if  the  ramp  slope  is 
not  large  enough  (condition  (ii)),  the  system  will  be  unstable. 

2.7  Comparison  with  Stability  Analysis  Using  State-Space  Averaging 

The  duty  cycle  changes  from  cycle  to  cycle.  This  variable  determines  when  the  switch  is  on  or 
off,  and  it  is  crucial  in  the  system  dynamics.  So  this  variable  needs  to  be  modeled  properly.  The 
averaging  method,  however,  treats  the  duty  cycle  as  a  continuous-time  variable. 

Next,  stability  analysis  using  the  averaging  method  is  briefly  summarized.  For  details,  the 
reader  is  referred  to  [15].  Let  the  duty  cycle  in  the  averaging  method  be  dc,  with  nominal  value 
Dc-  Also  let 


A 

-^ave 

:=  AonTIc  +  Aoff(1  - -Dc) 

(25) 

-Save 

:=  Rqn-Dc -1- -Boff(1  - -Dc) 

(26) 

E 

-^ave 

=  i?ON-Dc -1- -Doff(1  - -Dc) 

(27) 

-^ave 

~  ~^a.veEaveys 

(28) 

averaging  method  is  Aave,  a  constant  instead  of 

a  periodic  solution 

x^{t).  Linearized  around  the  nominal  operating  point  {x,u,dc)  =  (Aave,  (U,  U.)^  Uc),  the  system 
in  Fig.  3  has  the  following  linearized  continuous-time  dynamics 


X  ^  Aave®  +  -Bavefi  +  ((AqN  “  AoFF)A'ave  +  (RqN  “  Bo^-p)u)dc 

y  =  Cx  -|-  Du 

Vo  =  Es^eX  -I-  (F^ON  -  F'oFf)c?c 


(29) 

(30) 

(31) 


with 


dr  = 


y  _  y 

Vh-Vi  ~  h{d)T 

—y  —y 

Vh-Vi  ~  h{d)T 


(for  trailing-edge  modulation,  because  dc  =  Vh^Vi  ■ 
(for  leading-edge  modulation,  because  dc  =  VuZvi  ^ 


(32) 
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In  either  trailing-edge  modulation  {Si:  switch  on,  S2:  off)  or  leading-edge  modulation  {Si:  off, 


S2:  on).  Substituting  Eqs.  (29)  and  (32)  into  Eq.  (29)  gives 


X  ^  Ax  -1-  Bu  (33) 

where  (stage  Si  or  S2  can  be  either  on  or  off  stage) 

A  =  ^ave  +  ((^1  ^  ^2)^ave  +  {Bl  —  B2)u)  .  (34) 

h{d)T 

B  =  -Bave  +  ((^1  “  ^2)^ave  +  {Bl  —  B2)u)  ■  (35) 

h{d)T 

The  system  is  asymptotically  stable  if  all  of  the  eigenvalues  of  A  are  in  the  open  left  half  of  the 
complex  plane,  which  is  equivalent  to  all  of  the  eigenvalues  of  exp(^T)  being  inside  the  unit  circle. 
Erom  Eq.  (34),  exp(^T)  can  be  written  as  (for  either  trailing-edge  or  leading-edge  modulation) 


exp(.Ar)  =  e 


A2(T-d)+AT,d+ 


((■^1  — -^2  )-^ave  +  (-Bl  —B2)u)C 
h{d) 


(36) 


which  is  similar  in  form  to  $  in  Eq.  (16). 

Comparing  Eqs.  (36)  and  (16),  it  is  seen  that  the  averaging  method  agrees  with  the  sampled- 
data  method  under  the  following  conditions: 

1.  x^{d)  is  approximated  by  Xave  (true  if  the  periodic  solution  x^{t)  is  approximated  by  the 
equilibrium  Xave)- 

2.  Cx^{d-)  is  much  smaller  than  h{d). 

3.  T  is  very  small,  so  that  )-x°{d+))CT  ~  {x^{d~)  —  x^{d'^))CT. 

4.  The  matrices  Xi,  A2  and  {x^{d~)  —  x^{d~^))C  commute. 

Note  that  condition  4  is  unlikely  to  hold,  and  condition  2  generally  does  not  hold  for  current  mode 
control. 


2.8  Closed-Loop  Audio-Susceptibility  and  Output  Impedance 

Audio-Susceptibility  and  output  impedance  are  expressed  in  terms  of  transfer  functions  (frequency 
responses).  They  give  information  on  the  effect  of  (source  or  load)  disturbances  at  various  frequen- 
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cies  on  the  output  voltage. 


The  audio-susceptibility  is  derived  directly  from  the  linearized  sampled-data  model,  Eq.  (18). 
It  is 

Tos{z)  =  '^  =  E{zI (37) 

To  calculate  the  output  impedance,  add  a  fictitious  current  source  io  (as  perturbation)  in  parallel 
with  the  load.  Then  the  state  equations  in  Eqs.  (1)  and  (2)  are  replaced  by 


S'!  :  X 

=  Aix  +  BiiVs  +  BuVr  +  Bisio 

(38) 

S2  :  X 

=  A2X  +  B21VS  +  B22Vr  +  ^23*0 

(39) 

where  i?i3,i?23  G 

Since  io  is  used  as  perturbation,  the  nominal  value  of  io  is  0.  Similar  to  the  derivation  in 
Section  2.2,  the  new  linearized  sampled-data  dynamics  is 


Xn+l  ~  T  Tih^jj  -|-  T'2Vrn  T  r3io^ 

Von  —  EXji 


(40) 


where 


r  =  (  ^9  x-i  dg 

^  diem  ddn  dd„ 


{Xn,U,dn,ion)  =  {x^  {0),U,d,0) 


=  — x^{d+)^  t  da Bi^)  e^^‘^daB23  (41) 

Vo  Cx0{d-)-h{d)  Jo  Vo  ^  ^ 


and  ion  is  the  sampled  perturbed  output  current.  So  the  output  impedance  is 


Too{z)  =  ^  =  EizI  - 

to{z) 


(42) 


Given  a  transfer  function  in  z  domain,  say  T{z),  its  effective  frequency  response  [17,  p.  93]  is 
j'(giajT),  jg  valid  in  the  frequency  range  jenj  < 
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3  Discontinuous  Conduction  Mode  (DCM) 


3.1  Block  Diagram  Model 


There  are  three  stages  in  DCM  in  a  cycle.  The  first  two  stages  and  their  operation  are  the  same 
as  in  CCM.  The  system  is  switched  to  the  third  stage  when  the  inductor  current  reaches  zero. 
Within  the  third  stage,  =  Q.  A  block  diagram  model  for  the  PWM  converter  in  DCM  is  shown 
in  Fig.  9.  The  matrix  F  €  is  chosen  such  that  Fx  =  II-  The  remaining  notation  is  the  same 

as  in  Fig.  3. 


Figure  9:  A  block  diagram  model  for  PWM  converters  in  DCM 


3.2  Nonlinear  Sampled-Data  Model 

Consider  the  operation  of  the  PWM  converter  within  the  cycle,  t  €  [nT,  (n  +  l)r).  As  in  CCM, 
take  u  =  {vs,Vry  to  be  constant  within  the  cycle  and  denote  its  value  by  Un  =  {vsmVmY- 

Denote  by  nT  +  din  tke  switching  instant  when  y{t)  and  h{t)  intersect  within  the  cycle.  (The 
notation  din,  instead  of  di^n,  is  used  for  brevity.)  Denote  by  nT  +  d2n  the  switching  instant  when 
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the  inductor  current  reaches  zero.  The  three  stages  Si,  S2,  within  the  cycle  are  thus 


:  I  ^  for  t  G  [nT,  nT  +  din)  (43) 

I  Uo  —  r^ix 

<52  :  I  ^  ^  for  t  G  [nT  +  din,  nT  +  d2n)  (44) 

53  :  I  +  fortG[nr  +  d2n,(n  +  l)T)  (45) 

The  two  switching  conditions  (at  t  =  nT  +  din  and  nT  +  d2n)  are 

y{nT  +  din)  =  Cx{nT  +  din)  +  Dun  =  h{nT  +  di„)  (46) 

ZL(nT  +  d2n)  =  Fx{nT  +  d2n)  =  0  (47) 


As  in  CCM,  the  output  voltage  in  DCM  can  be  discontinuous  at  the  clock  time.  Again,  E  is 
used  to  denote  Ei,  E^,  or  {Ei  +  E^)/2  depending  on  which  value  of  output  voltage  is  of  interest. 
From  Eqs.  (43)-(47),  the  PWM  converter  in  DCM  has  the  following  sampled-data  dynamics: 


^n+1 


_  ^As{T  —  d2n)^^A2{d2n—d-\_n)^^Aldlnj.  _|_ 


eMdln-^)(^fjBiUn) 


pd2n  pT 

+  /  e^^^'^^^-^UaB2Un)  +  / 

J  d\ji  J  d2n 


Von  —  EXn 


g(Xn, Un, dn) 


Cx{nT  +  din)  +  Dun  —  h{nT  +  di„) 
Ex{nT  +  d2n) 


(48) 

(49) 


C'(e^i'^i"Xn  +  e^^^daBiUn)  +  Dun  -  h{din) 

J7(e^2(d2n-rfl„)(gAldl„3.^  ^  jdln  +  jd2n-dln  e^2^daB2Un) 


=  0 


(50) 


where  dn  =  (din,  d2n)^  Since  the  inductor  current  always  starts  from  0  at  the  beginning  of  a  cycle, 
another  explicit  constraint  is  Exn  =  ihu  =  0,  for  any  n.  So  the  dynamics  is  (iV-l)-dimensional 
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instead  of  iV-dimensional. 


3.3  Linearized  Sampled-Data  Dynamics 

Linearizing  Eqs.  (48)-(50)  and  using  the  notation  o  to  denote  evaluation  at  the  fixed  point  (x„,  Un,  dn) 
(x^(0),  u,  (di,  ^2)0)  has 


+  TUn 

EXri 


(51) 


where 


$  = 


r  = 


df 


dXn 

Of 


ddn 

dg 


ddr, 


dg 


dXn 

df 


dUn 

dg 


dUr 


df  dg 

V  n  J  / 


dXr,  ddr,  ddr,  dXr 


df  df  dg  1  dg 

'  1  a  j  ) 


du„  ddn  ddn  dUn 


gAs  (T— (i2 )  gA2  {d2  —di )  Q^idi 


^  gA3(T-d2) 


gA2(d2  x'^{d2  )  — 


CxO(df )  -  h{di)  0 

Fe^2(rf2-di)('^0(^-)  -  x0(4))  FxO(d^) 


C 

jp^A2(d2-di) 


oAidi 


gA3(T-d2)(gA2(d2-dl)  [ 

Jo  Jo 


Cd2—di 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

pT-d2 

e^^^daB2)+  /  e^^^daB^  (58) 


e^^^daBi  +  D 

p(^gA2(d2-di)  Jdi  ^^2-^1  gA2^^^52) 


(59) 


Since  the  inductor  current  ii,  starts  at  0  and  ends  at  0  in  any  cycle,  it  is  not  a  dynamic  variable. 
Thus,  the  dimension  of  the  dynamics  can  be  reduced  by  one.  Note  that  this  also  implies  that 
possesses  a  zero  eigenvalue,  i.e.  that  det[$]  =  0.  This  is  shown  in  detail  as  follows: 


det[d>]  =  det[e^3(T-d2)gA2(d2-di)('j_  r  e-^2(rf2-di)('^0(^-)  _ 
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(^)-l 


C 

p^A2{d2-d\) 


)e 


A\du 


=  det  [e^i  "*1  (T-da )  (^2-^1 )  ]  [ j 


C 

PQA2{d2-d\) 


x^{di)  —  x^{di)  e  ^ 


=  det[e^i'^ie"^3('^"'^2)gA2(<i2-rfi)]  . 


det  [I  — 


C {x^ {di  )  —  x^ {df ))  Ce  {x^ {d2  )  —  x^ {d^ )) 


Fe^2(rf2-rfi)(^0('^-)  _  a;0('^+)) 


Fx^{d7^ 


Fx^id^)  0 

-Fe^2(d2-<ii)(3,0('^-)  _  a;0(d+))  Cx^{d^)  -  h{di) 

(Cx^idi)  -  h{di))Fx^{d2) 


=  det [e^i (T-d2 ) gA2 id2-di ) j  j j  _ 


=  0 


here  *  signifies  an  irrelevant  term. 

The  stability,  audio-susceptibility  and  ontput  impedance  analysis  are  similar  to  the  case  in  CCM 
and  are  omitted. 


4  Illustrative  Examples 

Example  1  {Local  vs.  global  orbital  stability,  [9,  p.  90])  Consider  the  boost  converter  shown  in 
Fig.  10,  where  T  =  2g,s,  Vs  =  4F,  L  =  5.24/ii7,  C  =  0.2fiF,  R  =  160,  ki  =  —0.1,  k2  =  0.01, 
Vr  =  0.481^,  and  h{t)  =  ((j^)  mod  1). 

Althongh  the  control  scheme  is  neither  voltage  nor  current  mode  control,  the  circuit  can  still 
be  expressed  in  terms  of  the  block  diagram  model  in  Fig.  3,  with  state  x  =  (iL,vc)'- 


4i  = 

’  0  o' 

0  rU  _ 

4-2  = 

[0^1 
1  -1 
_  C  RC  _ 

Bi 

=  52  = 

1 

L 

0 

C  = 

—ki  —k2  D  = 

0  1 

El 

=  E2  = 

0  1 
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The  periodic  solution  calculated  using  Eq.  (8)  is  shown  in  Fig.  11.  The  closed- loop  poles 
calculated  from  Eq.  (16)  are  a($)  =  0.8  ±  0.45i,  which  are  inside  the  unit  circle.  So  the  periodic 
solution  is  locally  orbitally  stable.  The  magnitude  of  the  eigenvalues  is  0.9225,  indicating  that  the 
settling  time  to  the  steady  state  may  be  long.  For  example,  let  the  initial  state  be  {il-iVc)'  = 
(0.9,  8)'.  The  simulated  output  voltage  is  shown  in  Fig.  12.  After  the  transient,  the  state  trajectory 
goes  to  the  periodic  solution. 

The  averaging  method  also  predicts  local  stability  (but  not  local  orbital  stability).  From  Sec¬ 
tion  2.7,  the  closed-loop  poles  predicted  by  the  averaging  method  are 

a{A)  =  (-0.2759  ±  2.9276i)  x  10®  (60) 

and  =  0.7887  ±  0.5230i,  close  to  the  eigenvalues  found  above  using  by  the  sampled-data 

method. 

However,  the  circuit  is  not  globally  stable  and  is  described  in  [9]  as  being  unstable  based  on 
simulation.  Since  the  averaging  method  is  for  local  (small-signal)  analysis,  it  is  not  surprising  that 
it  cannot  predict  lack  of  global  stability. 


Figure  10:  System  diagram  for  Example  1 

Example  2  {Buck  converter  under  voltage  mode  phase-lead  control,  [15,  p.  346])  The  system 
diagram  is  shown  in  Fig.  13.  The  system  parameters  are  as  follows:  T  =  10/rs,  14  =  28F,  R  =  311, 
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L  =  bOfiH,  C  =  500//F,  Vr  =  5V,  Gco  =  3.7,  coz  =  10681  rad/ sec,  ujp  =  91106  rad/ sec,  and 
h{t)  =  4(|:  mod  1).  The  voltage  divider  gain  is  chosen  to  be  0.29465  (instead  of  1/3  in  [15])  to 
result  in  an  output  voltage  at  1517. 


Figure  13:  System  diagram  for  Example  2 


Let  the  state  x  =  {iL,vc,Xcy ,  where  Xc  is  the  state  of  the  error  amplifier.  In  terms  of  the 
representation  in  Fig.  3,  one  has 
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Solving  Eqs.  (9)  and  (13)  by  Newton’s  method  results  in  x°(0)  =  (4.3,15,-0.512)'  and  d  = 
5.36  X  10~®.  The  eigenvalues  of  the  closed-loop  system  can  be  obtained  from  (j[4>]  and  they  are 
0.8096  ±  0.1154  and  0.5973.  All  of  them  are  inside  the  unit  circle,  so  the  periodic  solution  x^{t)  is 
asymptotically  orbitally  stable. 

Example  3  {Boost  converter  under  current  mode  control  with  voltage  loop  closed,  [18])  The 
system  diagram  is  shown  in  Eig.  14,  where  fs  =  2bkHz,  Vg  =  28V,  R  =  11.211,  L  =  195/ri7, 
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C  =  2mF,  Ri  =  47.5kn,  R2  =  2.5ki},  Rs  =  0.81250,  Rf  =  72.2feO,  and  Cf  =  0.23nF.  The 
reference  voltage  Vr  and  thus  the  duty  cycle  are  varied.  The  system  is  analyzed  for  two  situations: 
without  slope  compensation  {h{t)  =  0),  and  with  slope  compensation  {h{t)  =  mod  1]). 

Let  the  state  x  =  {iL,vc,Vcfy ■  The  system  matrices  in  Fig.  3.1  are 
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First,  consider  the  system  without  slope  compensation.  The  duty  cycle  is  varied  (for  different 
Vr)  from  0.4  to  0.6  and  (t[$]  is  plotted  and  shown  in  Fig.  15.  One  eigenvalue  trajectory  crosses  the 
unit  circle  at  Dc  =  0.498.  The  other  two  remain  very  close  to  1.  For  Dc  >  0.498,  the  system  is 
unstable. 

Next,  consider  the  system  with  slope  compensation.  The  duty  cycle  is  varied  from  0.4  to  0.7 
and  cr[$]  is  calcnlated.  One  eigenvalue  trajectory  crosses  the  unit  circle  at  Dc  =  0.5845.  For 
Dc  >  0.5845,  the  system  is  nnstable. 

In  [18],  the  system  is  reported  to  be  unstable  for  Dc  >  0.454  without  slope  compensation  and 
for  Dc  >  0.61  with  slope  compensation. 

Example  4  (Boost  converter  under  current  mode  control,  [19])  Consider  the  same  boost  con¬ 
verter  as  in  Example  3,  with  the  following  differences:  here  fs  =  lt)QkHz,  Vr  =  1.817,  h(t)  =  0 
and  parasitic  resistances  in  the  switch,  diode,  inductor,  and  capacitor  are  modeled  and  are  given 
by  Rq  =  0.055n,  Rd  =  0.01111,  i?/  =  0.0311,  Rc  =  0.01211,  respectively.  With  these  choices,  the 
system  is  studied  in  [19]  and  shown  to  be  stable.  However,  both  state-space  average  model  and  an 
improved  model  [20]  predict  the  system  to  be  unstable. 

Let  the  state  x  =  (iL,vc,Vcfy ■  As  in  [19],  let  a  =  R/{R  +  Rc)  and  (3  =  aRc-  The  system 
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matrices  in  Fig.  3.1  are 
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The  closed-loop  poles  (c’'[$])  are  calculated  as  -0.3383,  0.9928  and  0.9994.  Thus,  the  sampled- 
data  approach  gives  a  better  indication  of  stability  than  the  averaged  model. 


5  Concluding  Remarks 


General  sampled-data  modeling  and  analysis  were  performed  for  closed-loop  PWM  DC-DC  convert¬ 
ers.  The  models  unify  past  work,  while  removing  simplifying  approximations.  They  apply  both  for 
current  mode  control  and  voltage  mode  control.  Linearized  models  were  derived  analytically,  and 
used  to  perform  associated  calculations,  among  which  are  general  results  on  stability.  Examples 
were  used  to  show  that  the  sampled-data  approach  will  reliably  predict  the  (local)  stability  of  a 
converter  for  which  averaging  or  simulation  predicts  instability.  Thus,  sampled  data  models  should 
be  viewed  as  an  important  tool  for  closed-loop  performance  validation.  It  is  hoped  that  this  work 
will  help  to  facilitate  further  applications  of  the  sampled-data  approach  in  power  electronics. 
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